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Abstract 

In previous work, the determination of uncertainty 
models via minimum norm model validation is based 
on a single set of input and output measurement 
data. Since uncertainty bounds at each frequency is 
directionally dependent for multivariable systems, this 
will lead to optimistic uncertainty levels. In addition, 
the design freedom in the uncertainty model has not 
been utilized to further reduce uncertainty levels. The 
above issues are addressed by formulating a min- 
max problem. An analytical solution to the min-max 
problem is given to within a generalized eigenvalue 
problem, thus avoiding a direct numerical approach. 
This result will lead to less conservative and more 
realistic uncertainty models for use in robust control. 

1 Introduction 

In applying multivariable robust control analysis and 
synthesis techniques to linear, time-invariant systems, 
as in for example [1], a set of plants as defined by a 
nominal and uncertainty model are required a priori. 
Nominal models are usually associated with a single 
“best” model, although what is considered “best” is 
debatable. Mathematical models derived from first 
principles are typically used as nominal models or 
sometimes identified from system identification exper- 
iments. In some cases where the physical conditions 
are not accurately or reliably known due to causes un- 
known or when simple models are desirable, it makes 
sense to require that the set of plants in question at 
least satisfy model validation conditions [2]- [6] with 
respect to available input and output measurement 
data. 

A recently proposed approach [4, 5, 6] to ob- 
taining uncertainty models from measurement data is 
based on searching for the smallest uncertainty set 
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that will validate all available data. The approach 
assumes a single set of input and output measurement 
data. For multivariable systems, this results in a 
particular input direction at each discrete frequency 
which in turn leads to optimistic uncertainty bounds. 

In addition, the need to admit model validating so- 
lutions implys a sufficient degree of freedom in un- 
certainty structure selection. This remaining freedom 
in a certain null space has not been utilized to reduce 
the minimum norm uncertainty bounds. In this paper, 
we investigate both issues. The problem formulation 
leads to a min-max problem and the conditions for a 
solution is given. 

We begin by first considering the standard PoA 
system shown in figure 1. Let the overall structured 
uncertainty be defined by the block diagonal matrices 

A = diag(A 1 , . . . , A r ); A® £ C™'*"', i = l,...,r(l) 

and the set of all block diagonal and stable, rational 
transfer function matrices be given by 

V = { A(-) G RHoo : A''(s 0 ) e C™' xn ',Vs 0 e C+} (2) 

where r and C+ denote the number of uncertainty 
blocks and the closed right-half plane, respectively [7]. 
The output error is given by 

e y := y Oy = y OP u (P, A)u (3) 

where the output, input, and upper linear fractional 
transformation of (P, A), are denoted by y, u, and 
tF u (P, A), and are variables in the z-domain so that 
they are complex vectors containing both gain and 
phase information. 

For model validation [4, 5, 6] with respect to input 
and output data u and y, the output error is set to 
zero so that from (3) 

y 0 P 22 U = T , 2iA(7 ^>PnA)~ 1 Pi 2 U (4) 

Observe from the left hand side of (4) that if the 
given nominal model P 22 and input cancels the mea- 
sured outputs (to within white noise), then A = 0 
validates the model. This is the case where a single 
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Figure 1: P<S>A structure in standard form. 

“optimal” model can be found such that the output 
errors are perfectly uncorrelated to the inputs. In 
this ideal case, a single “true” model is available and 
robustness issues with respect to model error do not 
exist. Otherwise, note that the nominal model error 
is directionally dependent on the inputs. Hence, with 
a different (directional) choice of inputs, the mininum 
norm model validating uncertainty computed will also 
be different. 

2 Problem Formulation 

To address the directional dependence of the mininum 
norm model validating uncertainty with respect to 
the test input, we assume that n u sets of inputs and 
corresponding outputs at each frequency are available 
and are arranged in the matrix forms 

U := [«i, . . . , «„,] ; Y := [t/i, . . . , t/„J (5) 

These sets of inputs and outputs can be formed from 
segments of a long identification data or independent 
identification experiments. Assuming a linear time- 
invariance system, a general input 

u = U a ( 6 ) 

will result in output 

y = Ya (7) 

so that the error equation for a general input signal is 

e y {a, A) = Ya <$f u {P, A)Ua (8) 

Note that a is a complex coefficient vector which 
help to span the n u dimensional vector space. The 
approach for handling the directional dependence is to 
search for the largest possible mininum norm model 
validating uncertainty over all possible test inputs. 

Consider model validating a linear combination 
of input/output data set, i.e., 

e y (a, A) = 0; Vo- (9) 

From ( 8 ), this condition is 


P 2 iA(I oPliAT 1 P 12 U a = (Y -&P 22 U)or, Va(10) 

Using the earlier approach and making the assump- 
tions as given in [4, 5], the A can be unscrambled to 


the form 
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where 



*( 0 -, <f>) 

:= Bee -(- P\\<f) 

( 12 ) 

y{<*A) 

:= Aa + <fr 

(13) 

and 



A = 

P+(Y^P 22 U) 

(14) 

B = 

P 12 U + P^P+iY oP 22 U) 

(15) 


4> £ M{P 2 i), where Af(-) denotes the null space. The 
variables depend on the arguments, a and <f>, and in 
the remaining development, these arguments will not 
be explicitly used for simplicity. 

To obtain lower bounds on the uncertainty com- 
ponents, partition in a conformal manner with respect 
to the uncertainty blocks the following: 

a- = col (x 1 , . . . , x T ) ; y = col (y 1 , . . . , y T ) (16) 
Due to the block diagonality of A, 

Ax = y, Ma O A l x l = y l , V(i, a) (17) 

and in addition 

ff’(A) = maxff’(A'), Vo- (18) 

i 

Consider the following Lemmas as given in [4, 5]. 

Lemma 1: 

For each component, x l 7 ^ 0, 

^(A f ) = ||A%>tt V(*». □ (19) 

Lemma 2: 

Let y { ' £ C m \ P £ C n \ P ± 0, 3A ?: £ C m ' xn ' : 

AV : = y\ a(A’ ) = |£|. □ 

The first, lemma gives a lower bound on the 
maximum singular value for the ith block diagonal un- 
certainty and the input/output data defined by a. In 
order to model validate the input/output relationship 
as defined by (17) for a given input vector function, 
and an output vector function, y l , the smallest norrned 
uncertainty is given by the lower bound in Lemma 
1. Of course if there are several input/output data, 
sets (corresponding to different, o-’s) to model validate, 
then the largest, lower bound with respect, to a is 
needed to validate all data, given. The second lemma. 





It can be shown that 


shows that the ith component uncertainty, A* , can 
be chosen to achieve this lower bound. Thus Lemma 
2 ensures that a minimum norrned uncertainty exists 
which will validate the given data. Again, if there are 
several input/output data sets the largest bound with 
respect to a is necessary. 

The problem is summarized in the following min- 
max problem: 


Problem 1 : max min J(a,<j>) 
« \4>e^(Pv2) 


1 1 :t/( q , ) 1 1 2 A 
IMt0)||2/ 


( 20 ) 


\\Aa + N4’\\l = fP T P(:=a(a,4’)> 0 
\\Ba + D'4’\\l = 4 T Q T Q4 '■= b(a , 4’) > 0 

win ly. the real column vector variable is 
' «r ] 

< = • l 

, V’i J 


(33) 

(34) 


(35) 


The physical interpretations for a and 4> are clear. The 
variable a represents all possible input disturbances 
in the system identification experiment whereas the 
variable 0 denotes the remaining freedom in the as- 
sumed uncertainty. Obviously, a larger number of un- 
certainty components in A will generally give greater 
freedom in 4> as realized by an increased dimension of 
the null space, Af. 

3 Solutions 

3.1 General Case 

In this section we consider the general min-max prob- 
lem defined in Problem 1. To avoid dealing with 
square roots, we consider the equivalent min-max of 
the ratio of norms squared. The problem considered 
is given as follows: 


max ( min 

a \ ip 

J(a,0(0)) 

04a + Jf'0 2 0 
• \\Ba + D4>\\ 2 J 

(21) 

where 




N H A = 

0 


(22) 

n h n = 

I 


(23) 

D = 

PuN 


(24) 


and 


P = 

[ A 

A ] = 

A r 
. A i 

04 i 
A r 

N r 

Ni 

N r 

(36) 

Q = 

[ Qi 

] = 

B r 

Bi 

B r 

D r 

Di 

D r 

(37) 


The positive definiteness of Q T Q is assumed from 
a physical argument. Observe that if b = 0 and 
a -jz 0, then J, the lower bound on the norm of model 
validating uncertainty is unbounded, which is to say 
a model validating uncertainty do not exist. Hence 
we assunnS that b > 0 V( ^ 0 so that Q T Q > 0 
or equivalently, Q has full column rank. The cost, 
function is then written as 


J(a,'4’) 


a(a, 0) 
v y ' > 0 
b(a,'4’) 


(38) 


As figure 2 illustrates, a solution to the min-max 
problem lies at a saddle point. The particular type 
of saddle point we seek must possess the following 
properties: it should be stationary with respect to £, 
local maximum with respect to a, and local minimum 
with respect to 0 . Specifically, 


The matrix N denotes an orthogonal basis for the null 
space Af and (22) shows its orthogonality with respect 
to the range space of the pseudo-inverse Pp. The free 
parameter 0 has the same dimension as the null space 
and is related to 0 by 


0 = 

N'4 ■ 



(25) 

Let 




A 

= 

-Ar + jAi 

^ (jrxn 

(26) 

N 

= 

N r + jNi 

^ £irxm 

(27) 

B 

= 

B r + jB{ 

^ (js xn 

(28) 

D 

= 

D r + jDi 

^ (js X m 

(29) 

a 

= 

a r T 0a ? ; 

e c mxl 

(30) 

4’ 

= 

Vv + j'4'i 

^ (jm X 1 

(31) 





(32) 



Figure 2: Geometry of a min-max solution: saddle 
point. 

dJ i-'./ dj 

de da T ’ d'4’ T 


0 stationarity (39) 



and 


d_ / dJ 
da \<9a T 
d_ ( dJ 
dip \dip T 


< 0 
> 0 


local maximum wrt a (40) 
local minimum wrt ip (41) 


By taking partial derivatives, it can be shown that the 
stationary condition reduce to 

[P t PoJQ t Q]£ = 0 (42) 


It can also be shown that the local maximum and 
minimum conditions reduce to the sign definiteness 
conditions 


PiPi&JQiQi < 0 (43) 

PfPi&JQZQi > 0 (44) 


Proposition 1: 

Let J k , !; k be the eigenpair of (42). Then the 
eigenpair is a solution to the min-max problem in 
(20) only if J k satisfy the inequalities (45). □ 


If more than one eigenvalue lie within the inequal- 
ity bound given in (45), this indicates that there is 
more than one saddle point and the smallest eigen- 
value should be chosen since the square root of these 
eigenvalues directly correspond to uncertainty lower 
bounds. From a computational standpoint, the real, 
symmetric generalized eigenvalue problem that needs 
to be solved is of dimension 2(n tl +dim(2V")), which will 
not likely be a problem for most practical applications. 

3.2 Special Case, <f> = 0 


Finally, it follows that the above negative and positive 
definiteness conditions in (43) and (44) respectively, is 
equivalent to the condition 


Consider the special case of <p = 0 or when Af(P 2 i) = 
{0}. The mini-max problem reduces to the single 
optimization problem: 


cr(PiL 1 T ) < J < g_(P2L 2 T ) 


(45) 


where Li are square, nonsingular matrices derived 
from the singular value decompositions of Qi [ 8 ], i.e, 

Qi = Ui £,■ Vf Li = ViZi, i=l,2. (46) 

Note that Q i, and Q 2 will have full column rank 
matrices since Q was earlier assumed to have full 
column rank on physical grounds. To recapitulate, 
the min-max problem defined by Problem 1 reduces 
to finding a and ip, such that J is stationary with 
respect to both a and ip, i.e., the equality (42) holds, 
and the curvature of J with respect to a is negative 
and is positive with respect to ip, i.e., the inequality 
(45) holds. 

To satisfy the conditions defined in (42) and (45), 
first note that both J and £ are unknowns in (42) so 
that this equation can be treated as a real, symmetric 
generalized eigenvalue problem. The functional de- 
pendence of both a and b and hence J on £ is not 
constraining as evident from premultiplying (42) by 
and solving for J . Hence a simple and direct solu- 
tion algorithm follows. First, find all stationary points 
by solving for all the eigenvalues and corresponding 
eigenvectors for the eigenvalue problem in (42) 


J k := 



e 


(47) 


where k = 1, . . . , 2 (n u + dim(7V")). Second, check to 
see if any of the Xj-’s lie within the inequality bound 
(45). The results are summarized in the following: 


Problem 2 : max J(a, 0) := 77 ^ 7 — ■ — yyp- 
« \\x(a, 0)|| 2 


I \Aa 1 1 2 
\\Ba\\ 2 


(48) 


Since A must be finite in satisfying the model vali- 
dation condition (9) for all nonzero a, B H B must be 
positive definite. With the nonsingular (and hence 
unique) coordinate transformation 


L t a = j3 

(49) 

the optimization problem reduces to 


MalC 

maxJ(a, 0 ) = max-— — — 
a a \\Ba 2 

(50) 

\\AL- t P3\\ 2 

— max - — - 

p mu 

(51) 

= a(AL~ T ) 

(52) 

The above development is summarized in 
ing. 

the follow- 

Proposition 2: 

The largest lower bound of the uncertainty with 
respect to all possible linear combinations of a set 
of n u input/output data is given as 

<t(A) > J(a,0) = a(AL T ). □ 

(53) 


Consider the special case where a specific lower 
bound is used instead of any linear combination. For 
example, let 6j denote the jth column of an n u by 
n u identify matrix. Then the jth input vector corre- 
sponds to a = 6j so that 



Uj =U6j; Vj = Y 6j (54) 

From (49) the corresponding [3 is 
Pi = L T 6j (55) 


so that from Proposition 2 


lk(«i.0)||2 


V 


11 / 3,112 


< a(AL T ) < cr(A); 


j = 1, • • • , n u 


(56) 


The expression on the left hand side of (56) denotes 
an uncertainty bound based on a single input/output 
data and is given in the earler works [4, 5]. 


3.3 Structured Uncertainty 

To compute norm bounds for unstructured uncer- 
tainty, the maximum singular value in (53) or the 
eigenvalues in (47) need to be computed for general 
and special case respectively. However, for structured 
uncertainty case, eigenvector f k and the maximizing 
principal vector [3 is needed to compute individual 
component uncertainty bounds for the general and 
special case respectively. 

For the special case, the norm bound on the ith 
component uncertainty is given from (19) as 

\\[AL~ t I3]% 


r(A*) > 


W 


(57) 


The superscript i in the right hand side of (57) de- 
notes the ith component of the column vectors. The 
principal vector, /?, comes from the singular value 
decomposition 

AL~ t = P • diag ( a(AL ~ T ), . . .) • [/?, . . ] T (58) 

Similarly, for the general case, the norm bound on the 
*th component uncertainty is also given from (19) as 

y(a k A(ip k ))h 


r(A*) > 


(59) 


where the subscript k in (59) denotes the smallest 
eigenpair solution that satisfies the min-max inequal- 
ity bounds. 


4 Conclusion 

The issue of directional dependence of multivariable 
uncertainty models is addressed along with the uncer- 
tainty design freedom beyond model validation. The 
problem is formulated as a search for a saddle point 
in a finite dimensional, complex space. Numerical 
and experimental validation and applications of the 
proposed solution is underway. It is expected that 
these results will lead to less conservative and more 
realistic uncertainty models for use in robust control. 
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